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ABSTRACT 

We consider observational tests for the nature of Ultraluminous X-ray sources (ULXs). These must distinguish 
between thermal-timescale mass transfer on to stellar-mass black holes leading to anisotropic X-ray emission, 
and accretion on to intermediate-mass black holes. We suggest that long-term transient behavior via the 
thermal-viscous disk instability could discriminate between these two possibilities for ULXs in regions of 
young stellar populations. Thermal-timescale mass transfer generally produces stable disks and persistent X- 
ray emission. In contrast, mass transfer from massive stars to black holes produces unstable disks and thus 
transient behavior, provided that the black hole mass exceeds some minimum value MBH.min- This minimum 
mass depends primarily on the donor mass and evolutionary state. We show that MBH.min ^ 50 M Q for a large 
fraction (> 90%) of the mass-transfer lifetime for the most likely donors in young clusters. Thus if long-term 
monitoring reveals a large transient fraction among ULXs in a young stellar population, these systems would 
be good candidates for intermediate-mass black holes in a statistical sense; information about the donor star is 
needed to make this identification secure in any individual case. A transient ULX population would imply a 
much larger population of quiescent systems of the same type. 
Subject headings: accretion, accretion disks — binaries: close — X-rays: binaries 



1. INTRODUCTION 

In the past few years high-angular-resolution observations 
with Chandra have revolutionized the study of X-ray bina- 
ries in nearby galaxies and have revealed whole populations 
of sources in a variety of g alaxy types (for a recent review 
see Fabbiano & White 2003). The detected X-ray fluxes have 
been combined with distance estimates to the host galaxies 
to infer the apparent X-ray luminosities of the sources, as- 
suming isotropic emission. The inferred X-ray luminosi- 
ties reveal a distinct class of sources: non-nuclear point 
sources with apparent X-ray luminosities above the Edding- 
ton limit for a ~ 10M Q black hole (> 2 x 10 39 ergs _1 ), of- 
ten referred to as ultraluminous X-ray sources (ULXs). The 
existence of suc h sources was f irst noted in EINSTEIN ob- 
servations (e.g., Fabbianol ll988l) . Sho rt-term variability de- 
tected in a number of them (see e.g., Fabbiano et al.1 l2003t 
iMatsumoto et alJ 120011) excludes the possibility of source 
confusion and strongly points towards accretion as the ori- 
gin of the X-rays. At present the majority of ULXs have 
been found mainly in young stellar populations and re- 
gions of recent star formation, althou gh a few have been 
identi fied in elliptical galaxies ( e.g., IColbert & Ptald l2002t 
Sarazin. Irwin. & Bregman 2001) with luminosities close to 
the lower end of the ULX range. 

If the apparent X-ray luminosities are indeed the true lu- 
minosities of the sources, their high values have very im- 
portant implications for their accreting compact objects. For 
sources with X-ray luminosities comparable or in excess of 
10 40 ergs _1 , the Eddington limit gives a lower limit on the 
mass intermediate between stellar (< 50 Mq) and supermas- 
sive (> 10 6 M Q ) black holes (BH). ULXs may thus suggest 
the existence of a new cla ss of compact objects: interm ediate- 
mass black holes (IMBH; Colbert & Mushotzky 1999). 

On the other hand it is still possible that the accreting 
co mpact objects in ULXs are of stella r mass (< 20 M ; 
see Be lczvnski. K alogera. & Bulikll2002l) . The high apparent 



X-ray luminosities can be explained in two different ways: 
(i) either the Eddington limit (rigorously derived for spher- 
ical accretion) is not relevant and i n fact can be exceeded 
(see Ruszkowski & Begelman 2003) or (ii) the apparent X- 
ray luminosities overestimate the tru e source luminosit ies be- 
cause the emission is anisotropic (IKing et alJ 12001). Al- 
though the theoretical basis for imposing the Eddington limit 
is somewhat unclear, there is strong support for it partly 
from observati ons of X-ray bursts from accreting neutron 
stars ( e.g.. iKuulkers et al.l2003ULewin. van P aradiis. & Taaml 
119951) and from the current unde rstanding of the evolut ionary 
history of wide binary pulsars (Webbink & Kalogera 1997) 
and Cygnus X-2 (where the compact object does not seem 
to have gained any significant amount of mass:lKing & RitterJ 
Il999t Kolb et al. 2000; Podsi adlowski & RappaportTTo OO). 

Anisotropic emission is probably associated with X-ray lu- 
minosities comparable to the Eddington limit. Binary sys- 
tems c an reach suc h high luminosities in two different situ- 
ations (King 2002): (i) thermal-timescale mass transfer typi- 
cally o ccurring when th e donor is more massive than the ac- 
cretor (IKing et alJl2001l). Cvgnus X-2 (IKing & RitterlfT99^) 
and SS433 (King, Taam, & Begelman 2000) may be exam- 
ples of this phase; (ii) X-ray transient outbursts, where the 
thermal disk instability governs the accretion behavior. The 
first possibility obviously requires donors more massive than 
black holes (> 3-5M Q ), and hence relatively young stellar 
environments, wh ereas the second must a pply to ULXs in old 
elliptical galaxies (Piro & Bildsten 2002). 

Although population studies suggest that a large 
fractio n of ULXs must be stel la r-mas s X-ray bina- 
ries ( Grimm. Gilfanov. & Sunvaevl Eool . some may 
contain IMB H. For the ULX i n M8 2, the very high peak 
luminosity ( Matsumoto et alJ 120011). t he q uasi-periodic 
oscillations ( Strohmaver & Mushotzkv 2003), and the 
detection of an isotropic nebula around it may point 
away from the anisotropic-emission possibility (although see 
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iKing & Poundsl2003l) . On the other hand, no ULXs are found 
inside den se clusters, where IMBH are expected not only to 
form jMiller & Hamiltonl2002l iPortegies Zwart & Mc Millan 
2000) but also remain, as they are much heavier than the 
average stellar mass in clusters (fast cluster disruption 
could help , but this is sue is beyond the scope of this paper; 
iGiirkan & Rasiol2003l) . Thus at present the physical origin of 
some of these sources is not clear, and there may be ULXs of 
both stellar and intermediate mass. 

In this Letter we suggest that long-term transient behavior 
(not just flux variability by factors o f a few to several) due 
to the thermal-viscous disk instability ( King. Kolb. & Burderi 
[l996t IKing &RitteH[T99l may distinguish the two possi- 
bilities for ULXs in regions of young stellar populations 
(< 10 8 yr). We show that one can define a minimum BH 
mass M BHi min f° r disk instability and thus transient behavior 
(§ 2). This minimum mass depends primarily on the mass 
and evolutionary stage of the donor star. We show that, for 
donor masses > 5 M (expected to be the most likely donors 
in young stellar environments), MBH.min ^ 50 M for a large 
fraction of the mass-transfer phase (> 90%). (§ 3). By con- 
trast, ther mal timesca le mass transfer is expected to be persis- 
tent iKing et all2001l) . Thus if long-term monitoring reveals 
a significant transient fraction among ULXs in a young stel- 
lar population, these systems would be good candidates for 
IMBH in a statistical sense; information about the donor star 
is needed to make this identification secure in any individual 
case. In § 4 we discuss the observational significance of this 
diagnostic and the connection to IMBH formation scenarios. 

2. MINIMUM BLACK HOLE MASS FOR TRANSIENT BEHAVIOR 

The thermal-viscous disk instability provides a currently 
accepted explanation for transient behavior in X-ray bina- 
ries. The instability causes the disk to undergo a limit cycle 
in which the central accretion rate passes through short high 
states (outbursts) and long low states (quiescence). This pic- 
ture was originally developed to explain dwarf novae but can 
be extended to soft X-ray transients by including the effects 
of disk irradiation (van Paradijs 1996; King, Kolb & Bur- 
deri 1996; Dubus et al. 1999; for reviews see Lasota 2001; 
Frank et al. 2002). The condition for transient behavior is 
that the disk surface temperature at its outer edge should lie 
below the hydrogen ionization temperature. This in turn re- 
quires the mean mass transfer rate to lie below a critical value 
M cr it (King. Kolb. & Burderi 1996). This depends primarily 
on the binary component masses Mbh, M2 and orbital period 
P, and to a lesser degree on the detailed vertical disk structure. 
The latter is of cou rse uncertain' in t his paper we use the form 
given in (eq. 32 in Dubus et a ll!999l) : 
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Althou gh the precise conditions assumed by Dub us etall 
( 1999) (in particular the central mass and vertical structure) 
probably cannot be extrapolated to all of the cases we shall 
consider, this equation gives an adequate idea of when tran- 
sien t behavior is likely. Using a somewhat simpler expres- 
sion King . Kolb. & Burderil {l996) first showed that the con- 
dition M < Merit translates into a minimum BH mass MBH.min 
required for the development of transient behavior. Similarly, 



equation (1) can be used to derive this minimum: 
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Our aim is to examine whether transient behavior favors a 
distinct BH mass range. We use mass transfer sequences cal- 
culated for a set of initial binary configurations (of varying 
orbital periods, black hole and donor masses) and we derive 
Merit for a given donor mass M2 and radius R2 (i.e., evolution- 
ary state). We then use the depen dence of Meri t on Mbh and 
disk radius given in (eq. 30 in Dubu s et all 1999ft and solve nu- 
merically for MBH.min by setting M cr it equal to the mass trans- 
fer rate found from our mass transfer sequences, for given 
(M2,/?2)- If the BH mass used in the mass transfer calcu- 
lations exceeds this minimum the system will be transient. 
(Note that for a given sequence, M depends most sensitively 
on the donor mass and evolutionary stage at the onset of mass 
transfer and not so much on the accretor mass.) We examine 
the range of values for MBH.min and whether transient behavior 
can be associated with certain types of BH accretors (§ 3.2). 

3. MASS TRANSFER SEQUENCES 
3.1. Stellar Evolution Code 

We calculate stellar models and mass-transfer sequences 
wit h an updated stellar evolution code described in detail 
in ittvanova et alJ l2003t iPods iadlowski. Rappap ort. & Pfahl 
2002). The current version has been modified to minimize 
numerical noise in the mass transfer calculations and ensure 
that the stellar and Roche lobe radii track one another during 
mass transfer. We use mixing length and overshooting param- 
eters of 2 and 0.25 pressure scale heights respectively. Since 
we are dealing with massive stars, we a ccount for mass los s 
due to stellar winds (rates adopted from lHurlev et all 2000). 
In calculating orbital changes we take account of both mass 
transfer and wind mass loss with the specific angular momen- 
tum of the mass-losing donor. We assume that any mass trans- 
fer above the Eddington rate is lost from the binary with the 
specific angular momentum of the accretor. 

We model mass transfer self-consistently following the 
donor response to the appropriate rate of mass loss. The 
mass-transfer rate M is calculated in an implicit manner, so 
that the donor radius R remains equal to the Roche lobe ra - 
dius 7?l (using Eggleton's approximation; Eggl etoijll983l) . 
We consider the radius-mass exponents of the Roche lobe 
(l = d lnR h /d InM and of the star itself ( = d InR/d InM in our 
solution method. The response of the Roche lobe to the mass 
transfer is solely a function of the mass ratio, whereas the re- 
sponse of the stellar radius depends on the mass transfer rate. 
For a given model, we tabulate values of C for a range of M 
values. We then identify the value of M for which the Roche 
lobe radius is equal to the stellar radius (predicted from the 
value of Q. In some cases, the solution for M is not unique; 
we then choose the lowest value to avoid large excursions in 
the rate. As the donor evolves, the stellar-radius response 
changes, so we recalculate the table of (,{M), if the predicted 
stellar radius differs from the calculated one by SlnR = 10~ 4 . 

3.2. Calculations and Results 
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To investigate the systematics of mass transfer and disk 
(in)stability we consider a large set of mass-transfer se- 
quences driven by Roche-lobe overflow: binaries with BH 
masses in the range 10- 1000M Q , donor masses in the range 
1 — 25 Mq. For each donor mass, we evolve single-star models 
to different evolutionary stages that cover most of the stellar 
lifetime: from the Zero- Age to the End of the Main Sequence 
(ZAMS and EMS), the Hertzsprung gap (HG), and through 
core helium burning. We consider the possibility of multiple 
mass-transfer episodes in the evolutionary history of each bi- 
nary and we evolve each of our models up to carbon ignition. 
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FIG. 1. — Two examples of mass-transfer sequences for a 1000 Mq BH 
with a 20 Mq donor. Mass transfer starts on the Zero- Age Main Sequence 
(solid) and the Base of the Giant Branch (dotted). The orbital period (top), 
mass transfer rate (middle), and derived minimum BH mass for transient be- 
havior (bottom) are shown as functions of time normalized to the total dura- 
tion of each mass transfer episode. 



Every evolutionary sequence of course starts and ends with 
transient behavior as the mass transfer rate rises from and re- 
turns to zero. For episodes of otherwise persistent mass trans- 
fer these transient windows form a very small fraction of the 
total mass transfer lifetime and have very low discovery prob- 
ability. In terms of the solutions for MBH.min, this means that 
at the start and end of every mass-transfer episode the mini- 
mum BH mass for transient behavior is very low and certainly 
enters the stellar-mass range. We eliminate these insignificant 
(due to their low discovery probability) transient epochs by 
excluding the first and last 5% of the mass-transfer lifetime. 

The behavior of two example mass-transfer sequences is 
shown in Figure 1 . These have been calculated for Mbh = 
1000 M Q and donors of 20 M Q at two evolutionary stages: 
unevolved (ZAMS) and at the base of the Giant Branch. Our 
results for the binary orbital period, mass transfer rate, and 
minimum BH mass for transient behavior are shown as a func- 
tion of time normalized to the total duration of each mass- 
transfer episode (~ 10 7 yr and ~ 2 x 10 3 yr, respectively). Ev- 
idently, for most of these episodes (90%-100% of their du- 
ration), transient behavior requires BHs in the intermediate- 
mass range MBH.min > 50 M© and not the stellar-mass range. 



In reality, for stellar-mass binaries the orbital separation is 
small enough that the radiation field of the 0,B donor is able 
to keep the disk ionized, and therefore stable. This effect is 
negligible for the much wider separations of IMBH binaries. 

For sequences with donors down to 10 M Q and 7M© re- 
sults are very similar to the 20 M Q sequences, with MBH.min > 
50 M© for more than 90% of their duration. We note that 
ULX luminosities would be reached in outbursts and would 
then probably reflect the Eddington luminosity rather than the 
mass-transfer rates shown in Figure 1 . 

For binaries with more evolved donors (during most of the 
short phase of core-helium burning when orbital periods are 
~ 10 yr), MBH.min values can be < 10M Q . In these cases ac- 
cretion disks are so large that their edges are cool and allow 
the disk instability to develop. However, such binaries are 
expected to be uncommon in young clusters for two reasons. 
First, they are too wide to survive stellar interactions; typ- 
ical interaction timescales are 10 5 yr for stellar densities of 
10 5 pc" 3 (typical of youn g stellar clusters in star-forming re- 
gions; see § 8.4 in Binnev & Tremaine 1987). Second, our 
numerical calculations show that their mass-transfer episodes 
are much shorter ( < 1 4 yr) than for MS donors (10 5 -10 7 yr). 
Therefore for the same formation probability, it is more diffi- 
cult to detect such short-lived X-ray phases. 

Sequences with 5M Q donors at different evolutionary 
stages show a qualitative change in behavior. They strad- 
dle along the dividing line between MBH.min values in the 
intermediate-mass and the stellar-mass range because of a 
close balance between two effects on MBH.min (eq- [2]): the 
orbital period and the mass trasfer rate at Roche-lobe over- 
flow. As a result a 5 M© donor at the ZAMS and beyond the 
MS gives MBH.min values in the stellar-mass, whereas the same 
donor filling its Roche-lobe close to the end of the MS leads 
to MBH.min values in the intermediate-mass range. Sequences 
with less massive donors always have MBH.min values in the 
stellar-mass range, i.e., < 50 Mq. 

Before drawing firm conclusions from our results, we ex- 
amine whether the derived MBH.min values are affected by the 
fact that we have used an IMBH for the calculations. We have 
repeated our sequences for a stellar-mass BH (10 Mq) and 
found no qualitative or significant quantitative differences, 
as is evident in Figure 2 for two example sequences with a 
10M Q donor. These deviations occur because the mass ra- 
tio is close to unity in one case, and lead to higher values of 
MBH.min for BH masses comparable to the donor masses. 

Given the qualitatively different results for donors more or 
less massive than ~ 5M Q , we consider their relative proba- 
bilities as Roche-lobe filling BH binary companions in young 
star-forming regions where most ULXs are found. From ba- 
sic results of stellar dynamics and preliminary calculations of 
our own (Ivanova, Kalogera, & Belczynski, in preparation) 
we find that relatively massive companions in are favored for 
a number of reasons: (i) BHs sink by dynamical friction to 
the center of young star-forming regions, as do massive stars, 
and therefore there is more of them in the BH's vicinity; (ii) 
massive stars have a higher cross section for capture by a BH 
and, if exchange into binaries is relevant, lower-mass objects 
are generally ejected in the interaction; (iii) such stellar inter- 
actions strongly favor orbital periods in excess of ~ 100 d, so 
even, if a low-mass companion were present at some point in 
the BH dynamical lifetime, it would not fill its Roche lobe in 
young clusters (requires orbital periods shorter than ~ 1 d). 

On the other hand, such relatively massive donors to stellar- 
mass BH could drive thermal-timescale mass transfer and 
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FIG. 2. — Comparison of orbital period (top), mass transfer rate (middle), 
and derived minimum BH mass for transient behavior (bottom), for mass- 
transfer sequences with 10 Mq donors at ZAMS for a 1000 Mq BH (solid) 
and a 10 Mq BH (open circles with reduced time resolution for best clarity 
of the figure). It is evident that results are insensitive to the BH mass used in 
the mass-transfer calculations. 



therefore produce persistent X-ray sources (King et al. 2001 ). 
Thus we conclude that, if a substantial fraction of ULXs prove 
to be transient, then IMBH accretors could be favored. 

4. DISCUSSION 

We have calculated mass transfer driven by relatively mas- 
sive stars (5-20M©) in BH binaries likely in young stellar 
environments, and derived a minimum BH mass for transient 
behavior that in the majority of relevant cases is in excess of 
50 Mq. This provides an observational diagnostic that could 
allow us to distinguish between stellar-mass (< 20 M Q ) and 
intermediate-mass BH binary models for ULXs. We note that 
in old populations of elliptica ls both classes of sources ar e 
expected to be transient (King 2002; Piro & Bildsten 2002). 



Hence transient behavior cannot be used as an observational 
diagnostic in old stellar systems (ages in excess of 10 s yr). 

So far there is only one candidate for a transient ULX. One 
(L x ^ 1-1 x 10 4 °ergs" ') is in a starburst galaxy NGC3628 
( Strickland et alJl2001l) . It may be associated with a ROSAT 
X-ray source (so the position is not well constrained) that 
faded below the sensitivity limit by a factor of more than 27 
and reappeared in Chandra observations. 

Current scenarios for IMBH involve formation in young 
stellar clus ters. One possibility in vokes repeated black hole 
mergers (Miller & Hamilton 2002), although gravitational ra- 
diation recoil ( Redmount & Rees 1989) could prevent this 
by ejecting merger products from the cluster. Another 
idea invokes runaway collisions of massive stars and even- 
tual collapse of the massive remnant (provided that stellar 
wi nds do not decrease the mass of the collision product; 
see lPortegies Zwart & McMillarJl200 0). Then an IMBH may 
form within the lifetime of the most massive stars (3 Myr), and 
it may acquire a binary companion within a cluster relaxation 
time after BH formation (at < 10 Myr), when stars as massive 
as ~ 20 Mq are still present. At that time binary separations 
would still be wide, favoring the formation of IMBH binaries 
with orbital periods longer than about 100 d. We are currently 
studying the dynamical evolution of an IMBH in young clus- 
ters and the characteristics of IMBH binaries and we expect 
to present our results in the near future (Ivanova, Kalogera, 
Belczynski 2003, in preparation). 

It is important to realize that the transient behavior we dis- 
cuss must occur on timescales far longer than an observer's 
lifetime. Therefore tests for transient behavior must be car- 
ried out in a statistical sense. Any quantitative statement ul- 
timately depends on the duty cycle and the outburst duration. 
It is also important to remember that such a detection would 
imply a much larger number ~ l/d> 100- 1000 of quiescent 
systems of the same type. 
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